Magnetic quantum oscillations in nanowires 
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Analytical expressions for the magnetization and the longitudinal conductivity of nanowires are de- 
rived in a magnetic field, B. We show that the interplay between size and magnetic field energy-level 
quantizations manifests itself through novel magnetic quantum oscillations in metallic nanowires. 
There are three characteristic frequencies of de Haas-van Alphen (dHvA) and Shubnikov-de Haas 
(SdH) oscillations, F = Fo/(l + 7) 3/2 , and F ± = 2F /\1 + 7 ± (1 + j) 1/2 \, in contrast with a 
single frequency Fq = Sphc/ (2ne) in simple bulk metals. The amplitude of oscillations is strongly 
enhanced in some "magic" magnetic fields. The wire cross-section area S can be measured using 
the oscillations as S = 4-7T 2 Sph 2 c 2 / ("ye 2 B 2 ) along with the Fermi surface cross-section area, Sf- 

PACS numbers: 72.15.Gd,75.75.+a, 73.63.Nm, 73.63.b 
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High magnetic fields have been widely used to explore 
the single particle spectrum of bulk metals. Historically, 
dHvA and SdH quantum oscillations in magnetic fields 
have provided an unambiguous signature and accurate 
quantitative information on the Fermi surface and the 
damping of quasiparticles Q] . Essential deviations from 
the conventional three-dimensional (3D) oscillations have 
been found in low-dimensional metals like 2D organic 
conductors 0> • At present conducting nanowires and 
nanotubes of almost any cross-section down to nanome- 
ter scale and of any length can be prepared with modern 
nano-technologies |J|. There are significant opportuni- 
ties for discovery of unique nanoscale phenomena arising 
from the dimension quantization. In particular, galvano- 
magnetic transport properties of nanowires have been the 
subject of many studies during last decades ||. Heremas 
et al. Q observed the semimetal-semiconductor phase 
transition in the magnetoresistance caused by the in- 
terplay between the electron cyclotron orbits, the size 
energy-level quantization and the inter-band transfer of 
carriers in Bi nanowires. Their magneto-conductance 
was theoretically addressed in the extreme ID limit 0. 
The Aharonov-Bohm-type oscillations of the magneto- 
conductance have been discovered in carbon nanotubes 
0, and connected with a metal-insulator transition 
caused by shifting of the van Hove singularities of the 
density of states [!j • More recently SdH oscillations were 
observed in arrays of 80 nm Bi-nanowires |10) and in 
200nm Bi-nanowires in first and second derivatives 
of resistance with respect to the magnetic field. There 
is a great demand for quantitative characterization of 
nanowires and analytical descriptions of the interplay be- 
tween dimension and field-induced energy-level quantiza- 
tions. 

In this Letter, we present the theory of magnetic quan- 
tum oscillations in long metallic nanowires in the longi- 
tudinal magnetic field, B, parallel to the direction of the 
wire z. We consider clean nanowires with the electron 
mean free path, I = vft comparable or larger than the 
cross size, R, but smaller than the nanowire length, L, 
which allows us to apply the conventional Boltzmann ki- 
netics. We also assume that the electron wavelength near 



the Fermi level is very small in the metallic nanowires, so 
that L ^S> I > R ^> 2nh/ (m*vp ), where vf is the Fermi 
velocity and m* is the band mass in the bulk metal. We 
find novel quantum oscillations of the magnetization and 
the conductivity caused by the interplay between mag- 
netic and dimension energy-level quantizations. 

Let us first calculate the magnetization M — —d£l/dB, 
where fi = -fcjjT^ ln[l + exp(-£ Q /(fc s T))] is the ther- 
modynamic potential, £ Q = E a — /1, E a is the single- 
particle energy spectrum and p is the chemical potential. 

Boundary conditions on the surface of the wire are 
not compatible with the symmetry of the vector poten- 
tial, A = B x r/2, so there are no simple analytical 
solution for E a in the magnetic field. However, one 
can overcome this difficulty in the quasi-classical limit, 
p fiw„, where uj s = ttvf/R, using the Tomonaga- 
like linearization of the energy spectrum 12] . Approx- 
imating the wire as an infinite round well one obtains 
E a = h 2 (k 2 m + k 2 )/(2m*). Here Hk is the continuous 
momentum along the wire, and discrete k nm are defined 
as zeros of the Bessel functions J\ m \(k nm R) = 0, where 
to = 0,±1,±2,... are the eigenvalues of z-component 
of the orbital momentum. In the quasi-classical limit 
J\ m \(k nm R) ~ {knmRy 1 / 2 cos(k nm R - 7t|to|/2 - tt/4), 
and k nm R = n(2n + l)/2 + n\m\/2 + ir/i with n = 
0,1,2,.... 

Hence, near the Fermi surface the spectrum is given 
by £, n ,m.k ~ huj s (2n + \m\ — Up) + h 2 k 2 / (2m*), which 
is identical to the spectrum in a parabolic "confine- 
ment" potential V(r) = m*uj 2 (x 2 + y 2 )/2 (here uf = 
fi/(Hu> s ) 3> 1). The major contribution to dHvA and 
SdH oscillations arises from the energy spectrum near 
the Fermi level, so we can replace the metallic nanowire 
with the confinement potential. In contrast with the 
original problem, the model Hamiltonian, H = (p — 
eA/c) 2 /(2m*) + V(r) + sfiBB has simple analytical eigen- 
functions, ^(r) oc exp(ikz)p^ exp(— p 2 /2)L^ (p 2 ) and 
eigenvalues 
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where p 2 = (m*u/ H) (x 2 + y 2 ) , uj 2 



ui c /2 1 lu c — eB/m*c, p B is the Bohr magneton, and a = 
{n, m, k, s} comprises all quantum numbers including the 
spin s = ±1. Using Eq.(l) and replacing negative m with 
— m — 1 one obtains 
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where pj = fi—tiu — (ibBs, (j, s = p — h(uj + uj + 
and e^ m (fc) = 2fi,rtw + hmoj ± + h 2 k 2 /(2m* 



k B T 

(2) 

\ — HbBs, 
). Sum- 



mations over n and m can be replaced by sums over 
r, r' = 0, ±1, ±2, ....±oo using twice the Poisson's formula 



and the variables x = 2um + uj^m and y 
in place of n and m, 
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We are interested in an oscillatory correction, £1 to 
the thermodynamic potential arising from the terms in 
Eq.(3) with nonzero r or r' . Introducing a new variable 
£ = x + h 2 k 2 /(2m*) — (J^, integrating by parts, extending 
the lower limit of £ down to — oo and taking routine inte- 
grals over k, J dkexp(iak 2 ) = (7r/ |a| ) - 1 / 2 exp[z7ra/(4|a|)] 

exp(y)] -1 = 



and over y 
-in/ sinh(7ra) 



£_/(k B T), fdyexp(iay)[l 
we finally obtain 
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where 



(4) 



A r (a:,y) 



B TL{2m*x/K) 1/2 cos[nrp B B/(hx)} /nry 
27rr 3 / 2 sinh[7r 2 fc s Tr/(fe)] COt V"^ 

(5) 

are oscillation amplitudes, and the summation formula 
J2 r ( z — r ) _1 = 7rcot(7rz) has been applied. 

Here and further we neglect quantum oscillations of 
the chemical potential. For the sake of transparency, we 
also neglect a damping of quantum levels by the impu- 
rity scattering in dHvA oscillations. We introduce this 
damping in the SdH effect (see below) neglecting quan- 
tum oscillations of the scattering rate 1/r. The quan- 
tum oscillations of p and 1/r could lead to a mixing 
of dHvA frequencies in multi-band metals as predicted 
and exp erimentally observed in several bulk compounds 
[l3L ll4Lll5l ITtH . However, they are negligible in the pres- 
ence of a field and size-independent "reservoir" of states 
(i.e. a sub-band with a heavy mass and the inter- 
band scattering. 

There are three characteristic frequencies 2lo and uj^ in 
the oscillating part of the magnetization A4 = —dCl/dB, 
rather then a single frequency u> c . 
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FIG. 1: Oscillating part of the magnetization versus the mag- 
netic field for relatively low fields and three temperatures. 
The resonance at u c — u> s /y/2 is due to a partial recovery of 
the energy-level degeneracy. 



The same frequencies are found in the conductance, a. 
The longitudinal conductivity is given by |T?H 
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where hw p = nk B T(2p + 1), p — 0,±1,±2,..., and 
^0 = 2nnk B T /h is the "frequency" of the "time"- 
dependent z-component of the vector potential, A z — 
iEcvQ 1 exp(—ivot), due to a longitudinal electric field 
E( ^ t < h/(k B T)). The static conductivity is cal- 
culated as the analytical continuation of this equation to 
v = iuq — > 0. The product of two GFs averaged over the 
random impurity distribution is factorized as the product 
of averaged GFs for a short-range scattering potential in 
absence of vertex corrections [01 , 



^ j(r)C(r^ 
^ ituSp - £, a 

n,m,k 



G s (r, r';w p ) 



where uJ p — lu p + lj p / (2\uj p \t). Then integrating the con- 
ductivity (r(r, vq) over the cross-section of the wire one 
obtains the conductance, 



a{v) = 



ie 2 k B T ^ 



h 2 k z /m* 



Lv m* (ifuJp - £ a )[ih(u P — vq) — £a\ 



itvJ p — £ Q 



(6) 



Integrating by parts the second diamagnetic term in 
Eq.(6) cancels the paramagnetic part at vq = 0. The rou- 
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tine analytical continuation [l9| of the remaining para- 
magnetic part yields the static conductance |20| in the 
limit v — » 0, 



0.2- 
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d^[SGf(0] 2 , (7) 



where G% (f) = [£-£„- ih/^r)}- 1 is the retarded GF 
and/® = [l + expC/^T)]- 1 . 

Summations over n and m are performed using twice 
the Poisson's formula, as in Eq.(3). The term with r = 
r' = yields the classical contribution, 



<>0 



Ee 2 T 



dxxx 



-x + (04 - x) 2 + ^ 2 /(4t 2 )) 1 / 2 ] 3 / 2 
((^ - x) 2 + ft 2 /(4r 2 )) 1 /2 



(8) 



after integrating over k and neglecting tempera- 
ture corrections. One can also neglect /i 2 /(4r 2 ) 
in the integral, Eq.(8) and obtain the conven- 
tional Drude conductance, <jq — Ne 2 r / (Lra*), 
where N = 8L{2m*) 1 / 2 fi 5 / 2 /(15h 3 nLj 2 ) is the to- 
tal number of electrons in the wire. Calculat- 
ing quantum corrections in a = ctq + a is sim- 
ilar to calculating of ft. Using the integrals 
/ dkk 2 exp(iak 2 ) = (i7r 1 /2/2)/|a|) 3 / 2 exp[i7ra/(4|a|)] and 
J dy exp(iay) cosh~ 2 (y) = — 7ra/ sinh(7ra/2) we obtain 

- \^V^ D / ±\ ( 7rr A 1 7rr(w+=F^ _ ) ^ 

a = y > B r (u,u )cos 

V ' ' \ to 2uj 4 



cos 

where 
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e 2 TkBTcos[irrfj,BB/(hx)] exp[— irr/ (2xr)] 



?i(2m*?ia;) 1 /2 r i/2 s inh[7r 2 fc B Tr/(fia;)] 

x cot — - . 

V 2x J 



(10) 



If the conventional dHvA frequency is high, Fq ^> B, 
three novel dHvA/SdH frequencies, F, .F of the wire 
can be estimated as F = B 2 /SB ~ ^B 2 \df / dB\/ (hj 2 ) 
with / = 2cj, w , 



and 



F-Fo/(l+ 7 ) 3 / 2 , 



F± =2F /|l + 7±(l+7) 



1/2 I 



(11) 



(12) 



where F = S F hc/ (2ne), 7 = 4w 2 /w 2 = 
47r 2 S , Fft 2 c 2 /(e 2 SB 2 ), S — ttR 2 is the cross-section area 
of the wire, and = n(m*) 2 v 2 F /fi 2 is the Fermi-surface 
cross-section area. They are related as F = F + F~\F + — 
F-\/(F+ + F-) 2 . 

Remarkably, both temperature and scattering damp- 
ing factors in Eqs.(5,10) depend on oj and oj rather than 
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FIG. 2: Oscillating magnetization for intermediate fields and 
three temperatures. Magic resonances are observed in many 
Fourier harmonics at low temperatures. 
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FIG. 3: Oscillating magnetization for high fields and two tem- 
peratures. 



on the cyclotron frequency uj c . Hence there are no con- 
straint on the value of the magnetic field imposed by 
those factors as soon as w s is large enough, u s > T, 1/r. 
In low fields, where 7 3> 1, all frequencies are much lower 
than F , F « F / 7 3 / 2 and F ± « 2F0/7. In high fields, 
where 7 « 1, two of them are about the same as Fq, 
F w F + ss Fq, while the third one appears to be much 
higher, F~ w 4F /7 ^> Fq. With respect to the Pauli 
paramagnetism and Landau diamagnetism in the bulk 
metal, amplitudes of quantum corrections in the magneti- 
zation and in the magnetic susceptibility, x, per unit vol- 
ume are about (fkj s / fi) 1 ^ 2 and \p/ '(hiJJs)] 1 / 2 , respectively, 
(to get these estimates we divide Ai by ttR 2 ). The rela- 
tive amplitude of quantum corrections in the conductance 
is (fkv s /[i) 5 / 2 , and about (huj s /fj,) 3 ^ 2 and (^//i) 1 / 2 in 
its first and second field derivatives, respectively. If we 
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take lo s about the same as u> c , the quantum corrections 
are much smaller than in the bulk metal, where they 
have the relative order of magnitude as [^JSpiOc)} 1 / 2 in 
M., [^/{fajJc)} 3 / 2 in x an d (huj c / fj,) 1 / 2 in a j|. However, 
there are some "magic" magnetic fields where the quan- 
tum corrections "explode". These are fields where the 
condition 2w/(w ± ) = (q + 2)/r is satisfied, so " cot" in 
Eqs.(5,f0) becomes infinite if q is an integer. In particu- 
lar, first harmonics with r = I become infinite if 

^ = q =111 (U) 

These magic resonances are clearly seen in Figs. 1, 2 at 
low temperatures, where we present numerical data for 
the oscillating part of the magnetization is 1000 

and we choose cos[7rr iibB/(Hx)] — 1). At high fields, 
lo c ^> uj s , the conventional dHvA pattern dominates, but 
the magic resonances are still there, Fig. 3. 

Let us elaborate more about the physical origin of the 
magic resonances. It is well known that the Landau lev- 
els are 5ei?/(27rcS)-fold degenerate in the bulk metal of 
the cross-section area S. The boundary conditions in 
the nanowire (approximated here by the confinement po- 
tential) remove the degeneracy, Eq.(l). Therefore the 
density of states at every level is reduced by a factor 
to 2 1 (ir 2 (MjJ c ), which explains the reduction of quantum 
amplitudes compared with the bulk metal. However, the 
magic resonance conditions partially restore the degener- 
acy of the spectrum, Eq.(l). For example, if uj c = ui s /-\/2, 
one obtains 2w = 3w s /y/2 and uj~ — uj s /\/2, so that 



E a = h 2 k 2 /(2m*)+huj s (6n+3\Tn\~m+3)/(2V2)+sfiBB, 
which is the same for all combinations of n and m with 
a fixed value of 6n + 3\m\ — m. Hence, compared with 
the amplitudes estimated above, the magic amplitudes 
are enhanced. The "anharmonic" corrections to the lin- 
earised energy spectrum in Eq.(l)imposed by the bound- 
ary conditions restrict their enhancement. 

It might be difficult to observe the novel oscillations 
in the magnetization of a single nanowire because its 
small volume, but they could be measured on bundles 
of nanowires. As far as SdH oscillations in nanowires 
[ToL ITU is concerned, their quantitative comparison with 
the present theory needs measurements in a wider field- 
range allowing for the reliable Fourier analysis. U sing the 
typical radius of J5i-nanowires R = 100 nm [(| llOL ITU 
and the Fermi surface cross-section area Sf = 10 13 cm -2 
[2ll | yields an estimate of hid s /kB ~ 50K with the car- 
rier mass m* — 0.1m e . Then the lowest temperature 
presented in Figs. 1, 2 is about 0.5 K with these param- 
eters. 

In conclusion, we have presented the theory of mag- 
netic quantum oscillations in clean metallic nanowires 
with simple Fermi-surfaces. We have found novel os- 
cillations caused by the interplay between size and 
field energy-level quantizations with three characteris- 
tic frequencies, calculated their amplitudes and identi- 
fied magic resonances, where the quantum corrections 
are strongly enhanced. Our findings suggest that one 
can measure both reciprocal and real space geometries of 
nanowires in a single measurement. 
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